A globally relevant stock of soil nitrogen in the Yedoma permafrost domain.

Nitrogen regulates multiple aspects of the permafrost climate feedback, including plant growth, organic matter decomposition, and the production of the potent greenhouse gas nitrous oxide. Despite its importance, current estimates of permafrost nitrogen are highly uncertain. Here, we compiled a dataset of >2000 samples to quantify nitrogen stocks in the Yedoma domain, a region with organic-rich permafrost that contains ~25% of all permafrost carbon. We estimate that the Yedoma domain contains 41.2 gigatons of nitrogen down to ~20 metre for the deepest unit, which increases the previous estimate for the entire permafrost zone by ~46%. Approximately 90% of this nitrogen (37 gigatons) is stored in permafrost and therefore currently immobile and frozen. Here, we show that of this amount, ¾ is stored >3 metre depth, but if partially mobilised by thaw, this large nitrogen pool could have continental-scale consequences for soil and aquatic biogeochemistry and global-scale consequences for the permafrost feedback.

Statistical upscaling of ecosystem CO2 fluxes across the terrestrial tundra and boreal domain: Regional patterns and uncertainties.

The regional variability in tundra and boreal carbon dioxide (CO2 ) fluxes can be high, complicating efforts to quantify sink-source patterns across the entire region. Statistical models are increasingly used to predict (i.e., upscale) CO2 fluxes across large spatial domains, but the reliability of different modeling techniques, each with different specifications and assumptions, has not been assessed in detail. Here, we compile eddy covariance and chamber measurements of annual and growing season CO2 fluxes of gross primary productivity (GPP), ecosystem respiration (ER), and net ecosystem exchange (NEE) during 1990-2015 from 148 terrestrial high-latitude (i.e., tundra and boreal) sites to analyze the spatial patterns and drivers of CO2 fluxes and test the accuracy and uncertainty of different statistical models. CO2 fluxes were upscaled at relatively high spatial resolution (1 km2 ) across the high-latitude region using five commonly used statistical models and their ensemble, that is, the median of all five models, using climatic, vegetation, and soil predictors. We found the performance of machine learning and ensemble predictions to outperform traditional regression methods. We also found the predictive performance of NEE-focused models to be low, relative to models predicting GPP and ER. Our data compilation and ensemble predictions showed that CO2 sink strength was larger in the boreal biome (observed and predicted average annual NEE -46 and -29 g C m-2  yr-1 , respectively) compared to tundra (average annual NEE +10 and -2 g C m-2  yr-1 ). This pattern was associated with large spatial variability, reflecting local heterogeneity in soil organic carbon stocks, climate, and vegetation productivity. The terrestrial ecosystem CO2 budget, estimated using the annual NEE ensemble prediction, suggests the high-latitude region was on average an annual CO2 sink during 1990-2015, although uncertainty remains high.

Spatial heterogeneity and environmental predictors of permafrost region soil organic carbon stocks.

Large stocks of soil organic carbon (SOC) have accumulated in the Northern Hemisphere permafrost region, but their current amounts and future fate remain uncertain. By analyzing dataset combining >2700 soil profiles with environmental variables in a geospatial framework, we generated spatially explicit estimates of permafrost-region SOC stocks, quantified spatial heterogeneity, and identified key environmental predictors. We estimated that Pg C are stored in the top 3 m of permafrost region soils. The greatest uncertainties occurred in circumpolar toe-slope positions and in flat areas of the Tibetan region. We found that soil wetness index and elevation are the dominant topographic controllers and surface air temperature (circumpolar region) and precipitation (Tibetan region) are significant climatic controllers of SOC stocks. Our results provide first high-resolution geospatial assessment of permafrost region SOC stocks and their relationships with environmental factors, which are crucial for modeling the response of permafrost affected soils to changing climate.

Large stocks of peatland carbon and nitrogen are vulnerable to permafrost thaw.

Northern peatlands have accumulated large stocks of organic carbon (C) and nitrogen (N), but their spatial distribution and vulnerability to climate warming remain uncertain. Here, we used machine-learning techniques with extensive peat core data (n > 7,000) to create observation-based maps of northern peatland C and N stocks, and to assess their response to warming and permafrost thaw. We estimate that northern peatlands cover 3.7 ± 0.5 million km2 and store 415 ± 150 Pg C and 10 ± 7 Pg N. Nearly half of the peatland area and peat C stocks are permafrost affected. Using modeled global warming stabilization scenarios (from 1.5 to 6 °C warming), we project that the current sink of atmospheric C (0.10 ± 0.02 Pg C⋅y-1) in northern peatlands will shift to a C source as 0.8 to 1.9 million km2 of permafrost-affected peatlands thaw. The projected thaw would cause peatland greenhouse gas emissions equal to ∼1% of anthropogenic radiative forcing in this century. The main forcing is from methane emissions (0.7 to 3 Pg cumulative CH4-C) with smaller carbon dioxide forcing (1 to 2 Pg CO2-C) and minor nitrous oxide losses. We project that initial CO2-C losses reverse after ∼200 y, as warming strengthens peatland C-sinks. We project substantial, but highly uncertain, additional losses of peat into fluvial systems of 10 to 30 Pg C and 0.4 to 0.9 Pg N. The combined gaseous and fluvial peatland C loss estimated here adds 30 to 50% onto previous estimates of permafrost-thaw C losses, with southern permafrost regions being the most vulnerable.

Widespread global peatland establishment and persistence over the last 130,000 y.

Glacial-interglacial variations in CO2 and methane in polar ice cores have been attributed, in part, to changes in global wetland extent, but the wetland distribution before the Last Glacial Maximum (LGM, 21 ka to 18 ka) remains virtually unknown. We present a study of global peatland extent and carbon (C) stocks through the last glacial cycle (130 ka to present) using a newly compiled database of 1,063 detailed stratigraphic records of peat deposits buried by mineral sediments, as well as a global peatland model. Quantitative agreement between modeling and observations shows extensive peat accumulation before the LGM in northern latitudes (>40°N), particularly during warmer periods including the last interglacial (130 ka to 116 ka, MIS 5e) and the interstadial (57 ka to 29 ka, MIS 3). During cooling periods of glacial advance and permafrost formation, the burial of northern peatlands by glaciers and mineral sediments decreased active peatland extent, thickness, and modeled C stocks by 70 to 90% from warmer times. Tropical peatland extent and C stocks show little temporal variation throughout the study period. While the increased burial of northern peats was correlated with cooling periods, the burial of tropical peat was predominately driven by changes in sea level and regional hydrology. Peat burial by mineral sediments represents a mechanism for long-term terrestrial C storage in the Earth system. These results show that northern peatlands accumulate significant C stocks during warmer times, indicating their potential for C sequestration during the warming Anthropocene.

Ecosystem carbon response of an Arctic peatland to simulated permafrost thaw.

Permafrost peatlands are biogeochemical hot spots in the Arctic as they store vast amounts of carbon. Permafrost thaw could release part of these long-term immobile carbon stocks as the greenhouse gases (GHGs) carbon dioxide (CO2 ) and methane (CH4 ) to the atmosphere, but how much, at which time-span and as which gaseous carbon species is still highly uncertain. Here we assess the effect of permafrost thaw on GHG dynamics under different moisture and vegetation scenarios in a permafrost peatland. A novel experimental approach using intact plant-soil systems (mesocosms) allowed us to simulate permafrost thaw under near-natural conditions. We monitored GHG flux dynamics via high-resolution flow-through gas measurements, combined with detailed monitoring of soil GHG concentration dynamics, yielding insights into GHG production and consumption potential of individual soil layers. Thawing the upper 10-15 cm of permafrost under dry conditions increased CO2 emissions to the atmosphere (without vegetation: 0.74 ± 0.49 vs. 0.84 ± 0.60 g CO2 -C m-2  day-1 ; with vegetation: 1.20 ± 0.50 vs. 1.32 ± 0.60 g CO2 -C m-2  day-1 , mean ± SD, pre- and post-thaw, respectively). Radiocarbon dating (14 C) of respired CO2 , supported by an independent curve-fitting approach, showed a clear contribution (9%-27%) of old carbon to this enhanced post-thaw CO2 flux. Elevated concentrations of CO2 , CH4 , and dissolved organic carbon at depth indicated not just pulse emissions during the thawing process, but sustained decomposition and GHG production from thawed permafrost. Oxidation of CH4 in the peat column, however, prevented CH4 release to the atmosphere. Importantly, we show here that, under dry conditions, peatlands strengthen the permafrost-carbon feedback by adding to the atmospheric CO2 burden post-thaw. However, as long as the water table remains low, our results reveal a strong CH4 sink capacity in these types of Arctic ecosystems pre- and post-thaw, with the potential to compensate part of the permafrost CO2 losses over longer timescales.

Tundra landscape heterogeneity, not interannual variability, controls the decadal regional carbon balance in the Western Russian Arctic.

Across the Arctic, the net ecosystem carbon (C) balance of tundra ecosystems is highly uncertain due to substantial temporal variability of C fluxes and to landscape heterogeneity. We modeled both carbon dioxide (CO2 ) and methane (CH4 ) fluxes for the dominant land cover types in a ~100-km2 sub-Arctic tundra region in northeast European Russia for the period of 2006-2015 using process-based biogeochemical models. Modeled net annual CO2 fluxes ranged from -300 g C m-2  year-1 [net uptake] in a willow fen to 3 g C m-2  year-1 [net source] in dry lichen tundra. Modeled annual CH4 emissions ranged from -0.2 to 22.3 g C m-2  year-1 at a peat plateau site and a willow fen site, respectively. Interannual variability over the decade was relatively small (20%-25%) in comparison with variability among the land cover types (150%). Using high-resolution land cover classification, the region was a net sink of atmospheric CO2 across most land cover types but a net source of CH4 to the atmosphere due to high emissions from permafrost-free fens. Using a lower resolution for land cover classification resulted in a 20%-65% underestimation of regional CH4 flux relative to high-resolution classification and smaller (10%) overestimation of regional CO2 uptake due to the underestimation of wetland area by 60%. The relative fraction of uplands versus wetlands was key to determining the net regional C balance at this and other Arctic tundra sites because wetlands were hot spots for C cycling in Arctic tundra ecosystems.

Nongrowing season methane emissions-a significant component of annual emissions across northern ecosystems.

Wetlands are the single largest natural source of atmospheric methane (CH4 ), a greenhouse gas, and occur extensively in the northern hemisphere. Large discrepancies remain between "bottom-up" and "top-down" estimates of northern CH4 emissions. To explore whether these discrepancies are due to poor representation of nongrowing season CH4 emissions, we synthesized nongrowing season and annual CH4 flux measurements from temperate, boreal, and tundra wetlands and uplands. Median nongrowing season wetland emissions ranged from 0.9 g/m2 in bogs to 5.2 g/m2 in marshes and were dependent on moisture, vegetation, and permafrost. Annual wetland emissions ranged from 0.9 g m-2  year-1 in tundra bogs to 78 g m-2  year-1 in temperate marshes. Uplands varied from CH4 sinks to CH4 sources with a median annual flux of 0.0 ± 0.2 g m-2  year-1 . The measured fraction of annual CH4 emissions during the nongrowing season (observed: 13% to 47%) was significantly larger than that was predicted by two process-based model ensembles, especially between 40° and 60°N (modeled: 4% to 17%). Constraining the model ensembles with the measured nongrowing fraction increased total nongrowing season and annual CH4 emissions. Using this constraint, the modeled nongrowing season wetland CH4 flux from >40° north was 6.1 ± 1.5 Tg/year, three times greater than the nongrowing season emissions of the unconstrained model ensemble. The annual wetland CH4 flux was 37 ± 7 Tg/year from the data-constrained model ensemble, 25% larger than the unconstrained ensemble. Considering nongrowing season processes is critical for accurately estimating CH4 emissions from high-latitude ecosystems, and necessary for constraining the role of wetland emissions in a warming climate.

Rapid carbon loss and slow recovery following permafrost thaw in boreal peatlands.

Permafrost peatlands store one-third of the total carbon (C) in the atmosphere and are increasingly vulnerable to thaw as high-latitude temperatures warm. Large uncertainties remain about C dynamics following permafrost thaw in boreal peatlands. We used a chronosequence approach to measure C stocks in forested permafrost plateaus (forest) and thawed permafrost bogs, ranging in thaw age from young (<10 years) to old (>100 years) from two interior Alaska chronosequences. Permafrost originally aggraded simultaneously with peat accumulation (syngenetic permafrost) at both sites. We found that upon thaw, C loss of the forest peat C is equivalent to ~30% of the initial forest C stock and is directly proportional to the prethaw C stocks. Our model results indicate that permafrost thaw turned these peatlands into net C sources to the atmosphere for a decade following thaw, after which post-thaw bog peat accumulation returned sites to net C sinks. It can take multiple centuries to millennia for a site to recover its prethaw C stocks; the amount of time needed for them to regain their prethaw C stocks is governed by the amount of C that accumulated prior to thaw. Consequently, these findings show that older peatlands will take longer to recover prethaw C stocks, whereas younger peatlands will exceed prethaw stocks in a matter of centuries. We conclude that the loss of sporadic and discontinuous permafrost by 2100 could result in a loss of up to 24 Pg of deep C from permafrost peatlands.

The positive net radiative greenhouse gas forcing of increasing methane emissions from a thawing boreal forest-wetland landscape.

At the southern margin of permafrost in North America, climate change causes widespread permafrost thaw. In boreal lowlands, thawing forested permafrost peat plateaus ('forest') lead to expansion of permafrost-free wetlands ('wetland'). Expanding wetland area with saturated and warmer organic soils is expected to increase landscape methane (CH4 ) emissions. Here, we quantify the thaw-induced increase in CH4 emissions for a boreal forest-wetland landscape in the southern Taiga Plains, Canada, and evaluate its impact on net radiative forcing relative to potential long-term net carbon dioxide (CO2 ) exchange. Using nested wetland and landscape eddy covariance net CH4 flux measurements in combination with flux footprint modeling, we find that landscape CH4 emissions increase with increasing wetland-to-forest ratio. Landscape CH4 emissions are most sensitive to this ratio during peak emission periods, when wetland soils are up to 10 °C warmer than forest soils. The cumulative growing season (May-October) wetland CH4 emission of ~13 g CH4  m-2 is the dominating contribution to the landscape CH4 emission of ~7 g CH4  m-2 . In contrast, forest contributions to landscape CH4 emissions appear to be negligible. The rapid wetland expansion of 0.26 ± 0.05% yr-1 in this region causes an estimated growing season increase of 0.034 ± 0.007 g CH4  m-2  yr-1 in landscape CH4 emissions. A long-term net CO2 uptake of >200 g CO2  m-2  yr-1 is required to offset the positive radiative forcing of increasing CH4 emissions until the end of the 21st century as indicated by an atmospheric CH4 and CO2 concentration model. However, long-term apparent carbon accumulation rates in similar boreal forest-wetland landscapes and eddy covariance landscape net CO2 flux measurements suggest a long-term net CO2 uptake between 49 and 157 g CO2  m-2  yr-1 . Thus, thaw-induced CH4 emission increases likely exert a positive net radiative greenhouse gas forcing through the 21st century.

A pan-Arctic synthesis of CH4 and CO2 production from anoxic soil incubations.

Permafrost thaw can alter the soil environment through changes in soil moisture, frequently resulting in soil saturation, a shift to anaerobic decomposition, and changes in the plant community. These changes, along with thawing of previously frozen organic material, can alter the form and magnitude of greenhouse gas production from permafrost ecosystems. We synthesized existing methane (CH4 ) and carbon dioxide (CO2 ) production measurements from anaerobic incubations of boreal and tundra soils from the geographic permafrost region to evaluate large-scale controls of anaerobic CO2 and CH4 production and compare the relative importance of landscape-level factors (e.g., vegetation type and landscape position), soil properties (e.g., pH, depth, and soil type), and soil environmental conditions (e.g., temperature and relative water table position). We found fivefold higher maximum CH4 production per gram soil carbon from organic soils than mineral soils. Maximum CH4 production from soils in the active layer (ground that thaws and refreezes annually) was nearly four times that of permafrost per gram soil carbon, and CH4 production per gram soil carbon was two times greater from sites without permafrost than sites with permafrost. Maximum CH4 and median anaerobic CO2 production decreased with depth, while CO2 :CH4 production increased with depth. Maximum CH4 production was highest in soils with herbaceous vegetation and soils that were either consistently or periodically inundated. This synthesis identifies the need to consider biome, landscape position, and vascular/moss vegetation types when modeling CH4 production in permafrost ecosystems and suggests the need for longer-term anaerobic incubations to fully capture CH4 dynamics. Our results demonstrate that as climate warms in arctic and boreal regions, rates of anaerobic CO2 and CH4 production will increase, not only as a result of increased temperature, but also from shifts in vegetation and increased ground saturation that will accompany permafrost thaw.